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Abstract. We investigate a nonstandard phase field model of Cahn-Hilliard type. The 
model, which was introduced in [16], describes two-species phase segregation and consists 
of a system of two highly nonlinearly coupled PDEs. It has been studied recently in [5] , [6] 
for the case of homogeneous Neumann boundary conditions. In this paper, we investigate 
the case that the boundary condition for one of the unknowns of the system is of third 
kind and nonhomogeneous. For the resulting system, we show well-posedness, and we 
study optimal boundary control problems. Existence of optimal controls is shown, and 
the first-order necessary optimality conditions are derived. Owing to the strong nonlinear 
couplings in the PDE system, standard arguments of optimal control theory do not apply 
directly, although the control constraints and the cost functional will be of standard type. 



1 Introduction 

Let f2 C IR 3 denote an open and bounded domain whose smooth boundary T has outward 
unit normal n, let T > be a given final time, and let Q := Q x (0, T) , E:=Tx (0, T) . 
In this paper, we study the following initial-boundary value problem: 

(e + 2p)p t + ^Pt - Ap = a.e. in Q, (1.1) 
5p t - Ap + f\p) = p a.e. in Q, (1.2) 

—— = 0, 7— = a(u — p) a.e. on E, (1.3) 
on an 

p(x, 0) = po(x) , p(x,0) = p (x) , for a. e. x G Q. (1.4) 

The PDE system fll.ip - fll.2l) constitutes a phase field model of Cahn-Hilliard type that 
describes phase segregation of two species (atoms and vacancies, say) on a lattice in the 
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presence of diffusion. It has been introduced recently in [16] and [5]; for the general physi- 
cal background, we refer the reader to [16]. The unknown variables are the order parameter 
p , interpreted as a volumetric density, and the chemical potential p . For physical reasons, 
we must have < p < 1 and p > almost everywhere in Q . The boundary (control) 
function u on the right-hand side of (II. 3p 2 plays the role of a microenergy source. More- 
over, e and 5 are positive constants, and the nonlinearity / is a double-well potential 
defined in (0, 1) , whose derivative /' is singular at the endpoints p = and p = 1 ; a typ- 
ical example is / = /1 + f?, with f 2 smooth and fi(p) — c(p log(p) + (1 — p) log(l —/?)), 
where c is a positive constant. 

The PDE system (jl.ip - fll.4p is singular, with highly nonlinear and nonstandard coupling. 
In particular, unpleasant nonlinear couplings involving time derivatives occur in (11. ip . 
and the expression f'(p) in ( II .2p may become singular. In the recent papers [5], [6], well- 
posedness and asymptotic behavior for t — > 00 and e \ of the system ( ll.ip -( fOj) were 
established for the case when the second boundary condition in (11. 3p is replaced by the 
homogeneous Neumann boundary condition dp/dn = 0; a distributed optimal control 
problem for this situation was analyzed in [JJ. We also refer to the papers [3] and jl], 
where the corresponding Allen-Cahn model was discussed. 

The paper is organized as follows: in Section 2, we state the general assumptions and 
prove the existence of a strong solution to the problem. Section 3 is concerned with the 
issues of uniqueness and stability. Section 4 then brings the study of a boundary control 
problem for the system (jl.ip - (jl.4p . We show existence of a solution to the optimal control 
problem and derive the first-order necessary optimality conditions, as usual given in terms 
of the adjoint system and a variational inequality. 

Throughout the paper, we make repeated use of Holder's inequality, of the elementary 
Young inequality 

ab < 7a 2 + — b 2 , for every a, b > and 7 > 0, (1.5) 
of the interpolation inequality 

y Q -y Q 

where p,g,r 6 [1,+ao], 6 G [0,1], and - = h , (1.6) 

r p q 

and, since dimf2 < 3, of the continuity of the embeddings H 1 ^) C L q (Q) for 1 < q < 6, 
where, with constants C q > depending only on Q , 



1 1,9(0) 



<C q \\v\\ mm \/veH\n), l<q<Q, (1.7) 



and where the embeddings are compact for 1 < q < 6 . We also use the Sobolev spaces 
H S (Q) of real order s > and recall the compact embeddings H S (Q) C i/ 1 (fi) and 
H s (fl) C C(Q) for s > 1 and s > 3/2, respectively, and, e.g., the estimate, with a 
constant > depending only on Q , 



v 
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2 Problem statement and existence 

Consider the initial-boundary value problem fll.ip - fll.4p . For convenience, we introduce 
the abbreviated notation 

H = L 2 (n), V = H 1 ^), W = {we H 2 (n) : dw/dn = on T} . 

We endow these spaces with their standard norms, for which we use self-explaining nota- 
tion like || ■ \\v ; for simplicity, we also write \\-\\h for the norm in the space H x H x H . 
Recall that the embeddings W C V C H are compact. Moreover, since V is dense 
in H, we can identify H with a subspace of V* in the usual way, i.e., by setting 
(u,v)v*,v = (u,v)h for all u G H and v G V, where (• , -)v*,v denotes the duality 
pairing between V* and V . Then also the embedding H C V* is compact. 

We make the following assumptions on the data: 

(Al) / = A + / 2 , where A e C 2 (0, 1) is convex, A e C 2 [0, 1] , and 

limA(r) = -cx), limA(r) = +oo. (2.1) 

r\0 r/*l 

(A2) Po eW, f'( Po )eH, fi EV,and 

< p (x) < 1 VxGO, /i > a.e. in fi. (2.2) 

(A3) n G ^(O.T;!, 2 ^)), and u > a.e. on E. 

(A4) a G L°°(r) , and a(x) > ao > for almost every x G T. 

Notice that (A2) implies that po £ C(f2) and, thanks to the convexity of A > a l so that 

f(po)&H. 

The following existence result resembles that of Theorem 2.1 in [5]. 

Theorem 2.1 Suppose that the hypotheses (Al)-(A4) are satisfied. Then the system 
(jl.lj) — (jl.4p /ios a solution (p, /i) swc/i i/iai 

p G ^°°(0, T; if) n ^(0, T; V) PI L°°(0, T; VT), (2.3) 

p G H\0 : T;H) nC°([0,T\;V) nL 2 (0,T;F 3/2 (fi)), (2.4) 

A(p)GL°°(0,T;iJ), (2.5) 

< p < 1 a.e. in Q, p > a.e. in Q. (2.6) 

Remark 2.2 The # 3 / 2 space regularity for p is optimal due to the L space regularity 
of u given by (A3). Nevertheless, both equation f 1 1.11) and the boundary condition for 
p contained in (j 1.3 ft can be understood a.e. in Q and a.e. on S , respectively, and the 
standard integration by parts is correct, as we briefly explain (so that we can both refer to 
that formulation and use integration by parts). In principle, one can replace the equation 
and the boundary condition by the usual variational formulation, namely 

/ [(e + 2 p)pi + ppi] v alx + / Vp-Vvalx+ I a(p — u)v da — 
Jn Jn Jr 
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(where da stands for the surface measure) for every v G V , a.e. in (0, T) , or an integrated- 
in-time version of it. This implies that (II. ip is satisfied in the sense of distributions, 
whence Ap belongs to L 2 {Q) by comparison, and the equation can be understood a.e. 
in Q , a posteriori. The last regularity (12 .4p of p and the condition A/x G L 2 (0,T; 7 2 (ft)) 
just observed also ensure that the trace has a meaning in the space L 2 (0,T; 7 2 (t)) 
due to the trace theorem [T5l Thm. 7.3] (we just observe that the space E _1//2 (ft) that 
enters such a result is larger than L 2 (ft)), so that the boundary condition can be read 
a.e. on £ . 

Proof of Theorem 2.1. The proof follows closely the lines of the proof of Theorem 2.1 in 
[5] , where a homogeneous Neumann boundary condition for /x was investigated. 

Step 1: Approximation. We employ an approximation scheme based on a time delay in 
the right-hand side of (II. 2p . To this end, we introduce for r > the translation operator 
% ■ ^(0, T; H) -> L\0, T; H) , which for v G L\0, T; H) and almost every t G (0, T) is 
defined by 

(T T )(t) := v(t - r) if t > r, and (T T ){t) := /x if t < r . (2.7) 

Now, let iV G IN be arbitrary, and r := T/N. We seek functions (p r ,p r ) satisfying 
(I2.3l) - (l2.6p (with (p, p) replaced by (p T ,p T )), which solve the system 

(e + 2 p T )/x[ + /x r p[ - A/x r = a. e. in Q, (2.8) 

Spl - Ap T + f'(p T ) = T T (i T a. e. in Q, (2.9) 
0p T dp7 

— — = , — — = a(u — a T ) a. e. on E, (2.10) 
on on 

p T (x, 0) = po(x) , /x r (x, 0) = po(x) , for a. e. ifQ. (2-H) 

We note that Remark 2.2 also applies to the approximating problem. To prove the 
existence of a solution, we put t n := nr, I n :— [0, t n ], 1 < n < N, and consider for 
1 < n < N the problem 

(e + 2p n )/x™ + p n p™ - Ap" = a. e. in ft x I n , (2.12) 

p n (0) = p a.e. in ft, = a(u - p n ) a. e. on T x /„, (2.13) 

5p™ - Ap n + /'(p n ) = r r p n_1 a. e. in ft x J n , (2.14) 

p n (0)=p a.e. in ft, ^- = 0, a.e. on T x J n . (2.15) 

on 

Notice that the operator T T acts on functions that are not defined on the entire interval 
(0, T) . However, its meaning is still given by ( 12 .7p if n > 1 , and for n = 1 we simply put 

TrpT 1 - 1 = p • 

Clearly, we have (p r , p r ) = (p N , p N ) if (p^, p w ) exists. We claim that the systems (12.121) - 
(I2.15P can be uniquely solved by induction for n = 1, N , where, for 1 < n < N , 

p n e W l,ao( In . H) n H l( In . y) n L oo (4 . w)> (2 16) 

p n G ^(4; £T) n C°(7 re ; y) n L 2 (I n , if 3/2 (ft)), (2.17) 
< p n < 1 a. e. in ft x I n , p™ > a. e. in ft x 7 n . (2.18) 
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To prove the claim, suppose that for some n G {1,...,A} the problem (I2.12p ~ p.15p 
has a unique solution satisfying (12.160 — (12.1 8p . where the index n is replaced by n — 1 . 
Then it follows with exactly the same argument as in the proof of Theorem 2.1 in [5] that 
the initial-boundary value problem (I2.14p . ( I2.15P has a unique solution p n that satisfies 
(I2.16P and the first inequality in (I2.18p . Substituting p n in (I2.12p . we infer that the linear 
initial-boundary value problem (I2.12p . (I2.13P has a unique solution p n satisfying ( I2.17p . 
Notice here that the regularity of p n follows from the fact that u G if 1 (0,T; L 2 {T)) . 
It remains to show that p n is nonnegative almost everywhere. To this end, we test (I2.12p 
by — (/i n )~, where (p n )~ denotes the negative part of p n . Using integration by parts 
and the boundary condition in ( I2.13p . we obtain the identity 

\ I I ^-({e + 2p n ) \{fi n y\ 2 )dxds+ {I \V{p n )~\ 2 dxds 
- Jo Jn dt V / J J n 



+ / / a | (/i n ) | da ds + / au (/i n ) da ds = . 
Jo Jr Jo Jr 

From the fact that p n , p , p , a , u are all nonnegative, we infer that 

\{p n )-{t)\ 2 dx < [ (e + 2p n (t)) \(p n )-(t)\ 2 dx 



< (s + 2p ) | Mo \ 2 dx = . 
Jn 

Hence, {p n )~ = 0, i.e., p n > a. e. in Q x I n , and the claim is proved. 

Step 2: A priori estimates. Now that the well-posedness of the problem 
(I2.8p - (l2.1ip is established, we perform a number of a priori estimates for its solution. 
For the sake of a better readability, we will omit the index r in the calculations. In 
what follows, we denote by C > positive constants that may depend on the data of the 
system but not on r. The meaning of C may change from line to line and even in the 
same chain of inequalities. 

First estimate. Since d t ((e/2)p 2 + p/i 2 ) = ((e + 2p)p t + pp t ) p, testing of (12.81) by p 
yields, for every t G [0, T] , 

/ ^ p 2 + pp 2 j (t) dx + / / | V p\ 2 dx ds + / I a p 2 da ds 

= J ^- /ig + poPo^jit) dx + y J a up dads, 

whence, using Young's inequality and (A2)-(A4), we can conclude that 

||M||L oo (0,T;//)nL 2 (0,T;y) < C • (2-19) 



Second estimate. Next, we test (I2.9P by p t . Applying (I2.19p . recalling the fact that 
f(po) G H , and invoking Young's inequality, we easily see that 
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\p\\h 1 {0,T;H)C\L°°{0,T;V) + || f(p) \\ L°°(0,T;IA(p)) < C . 



(2.20) 



Third estimate. We rewrite Eq. (12. 9p in the form 

~ A P + f[(p) = -Sp t - f' 2 {p) + %p 

and observe that the right-hand side is bounded in L 2 (Q) . Hence, applying a standard 
procedure (e.g., testing by f[(p)), and invoking elliptic regularity, we find that 



\p\\l*{o,T;W) + llii(p)IU 2 (Q) < C. 



(2.21] 



Fourth estimate. We differentiate Eq. (12. 9p formally with respect to t and test the 
resulting equation with p t (this argument can be made rigorous, see [5]). Since, owing to 
the convexity of f\ , f"(p) is nonnegative almost everywhere, we find the estimate 



~IIM*)I| 2 



\Vp t \ 2 dxds < - ||Ap - f'i{po) +A*o||h 



+ max I/2 (p) | / / \p t \ 2 dxds + / / (d t T T p) p t dxds 
o<p<i " n J Jn JoJn 



< C 



t-T 



Jn 



Pt{ s ) Pt{s + t) dx ds . 



(2.22) 



In order to estimate the last integral, we substitute for fi t , using Eq. (12. 8p . It follows, 
using integration by parts: 



t-T 



Pt Pt(- + T ) dx ds 



Jn 

t-T 



t-T 



Jn L 



t-T 



V/i 

e + 2p 

1 



e + 2p 



PtPPt{- + r 



Jn^ + Zp 

2p t {- + t 



(e + 2pf 
dx ds 



(Afi — p p t ) pt(- + t) dx ds 
V/i • Vp 



e + 2p 



(u — p) pt(- + t) da ds . 



(2.23) 



Exactly as in the proof of Theorem 2.1 in [S], the domain integral in the second and third 
lines of (I2.23P can be estimated from above by an expression of the form 

\JJjVp t \ 2 dxds + C(l + j\\p{s)\\ 2 v \\p t {s)\\ 2 H dxds). (2.24) 

Observe that, owing to the inequality (I2.19p . the mapping s \-> \\p(s)\\y belongs to 
^(0,T). 
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Finally, we estimate the boundary term in the last line of Eq. ( 12. 23p . To this end, 
recall that by the trace theorem there is a constant cn > , independent of r , such that 
IMU 2 (r) < cq \\ v \\v for all v G V . Moreover, we have p > and a G L°°(r) . Therefore, 
we obtain that 

*" r ' a 

(u — ji) p t {- + r) da ds 



o Jv 



e + 2p 



t-T 



< C / \\p t (s + T)\\ L 2 (r) (\\u(s)\\ L 2 {r) + \\p(s)\\ L 2 (r) ) ds 







t-T 



< C / ||ft(a + T)|| v (lk( S )|U 2( r) + ||A*0Ollv) ds 

Jo 

< \f\\pt(s)\\ 2 v ds + C. (2.25) 

Now we may combine the estimates fl2.22p - fl2.25p and employ Gronwall's inequality to 
conclude that 

\\pt\\L°°(0,T;H)nL 2 (0,T;V) < C . (2.26) 

The same argument as in the derivation of (12.221) then shows that also 

lbt|U°°(0,T;W) + \\f[(p)\\L°°(0,T;H) < C. (2.27) 

Fifth estimate. We test equation (12. 8 j) by /i t . Formal integration by parts (this can be 
made rigorous), using (A3), the trace theorem and Young's inequality, yields: 

< C + / aup t da + / \fi p t p t \ dx ds 

Jo Jv Jo Jq 

< C + au(t)p(t)da — / au t p,da ds + / \pp t fit\ dxds 

Jv Jo Jv Jo Jq 

t rt 



< 1- 7lK*)||y + / !IM s )llv ds + / / \fi\\p t \\p t \ dx ds 

7 Jo Jo Jn 

< - + l\Ht)\\v + / Hs)||^d5 +1/11^)111,^ 
7 Jo 2 Jo 



/ 



+ CI \\pt(s)f L A m \\p{s)\\\ im ds 
Jo 

< ^- + 1 \\^t)\\ 2 v +if \\^( s )\\ 2 H ds 

7 ^Jo 

+ C[\l + \\p t {s)f v ) \\p(s)r v ds. (2.28) 
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Hence, using (I2.26p . choosing 7 > sufficiently small, and invoking Gronwall's lemma, 
we can conclude that 

Wf J >\\H 1 (0,T;H)nL°°(0,T;V) < C . (2.29) 

Sixth estimate. Since < p < 1 a. e. in Q, and using (I2.26p . (I2.29P and the continuity 
of the embedding V C L 4 (f2) , we can estimate as follows: 

|| (e + 2p) fl t + fip t || L 2(Q) < C \\p t \\ L 2(Q) + ||/i||z,°o( ,T;L 4 (C)) ||pJl 2 (0,T;L 4 (C)) 

< C (ll^t||L 2 (Q) + |HU°°(0,T;V) ||Pt||L 2 (0,T ; y)) < C . (2.30) 

Comparison in ( 12 .8p then shows the boundedness of Ap in L 2 (Q) , and it follows from 
(12. 8p . (A3) and standard elliptic estimates that also 

IHIl 2 (0,T;#3/2( Q )) < C . (2.31) 

Step 3: Conclusion of the proof. Collecting all the above estimates, it turns out that 
there is some sequence Tk \ such that 

p rfc — > p weakly star in 

H\0, T; H) n L°°(0, T; 1/) n L 2 (0, T; H 3/2 (n)) , 
p Tfc -»■ p weakly star in W 1 ' 00 ^, T; H) n if^O, T; V) n L°°(0, T; W) , 
/iV"') -> i weakly star in L°°(0, T; if) . 

Thanks to the Aubin-Lions lemma (cf., [T4J Thm. 5.1, p. 58]) and similar results to be 
found in [T71 Sect. 8, Cor. 4], we also deduce (recall that even ii 3//2 (f2) is compactly 
embedded into V) the strong convergences 

p Tfc -> p strongly in C°([0, T]; if) n L 2 (0, T; 1/) , 
p Tfc -»■ p strongly in C°([0, T]; ^) 

and the Cauchy conditions (II. 4p as a consequence. In particular, employing a standard 
monotonicity argument (cf., e.g., [H Lemma 1.3, p. 42]), we conclude that < p < 
1 and £ — f[ (p) a. e. in Q . The strong convergence shown above also entails that 
f2i.P Tk ) ~^ f^ip) strongly in C°([0, T]; H) (because f 2 is Lipschitz continuous), and that 
T Tk \x Tk —> p strongly in L 2 (Q). 

Now notice that the above convergences imply, in particular, that 

p Tk -> p strongly in C°([0, T]; L 6 (fi)) , 

p[ fc -> pt weakly in L 2 (0, T; L 4 (fi)), 

p Tfc -»■ p strongly in L 2 (0, T; L\fl)) , 

p[ fc — » P( weakly in L 2 {Q) . 
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From this, it is easily verified that 

p Tk p Tk — >■ p, pt weakly in L 1 (0, T; H), 

p Tk p\ k ->■ pp t weakly in L 2 (0, T; L 3/2 (fi)). 

Now, we are ready to take the limit as k — > oo in fl2.8l) - fl2.10l) (written for r = r fe ). 
Precisely, we can do that as far as p is concerned, while it is easier to take the limit in 
the variational formulation of (12.81) that accounts for the boundary condition (the same 
as mentioned in Remark 2.2), or in the following integrated-in-time version of it 

f I [(e + 2p T )p T t + p T p T t ] vdxdt+ [ [ V p T -Vvdxdt 
Jo Jn Jo Jn 

a(p T - u)vdadt = for every v G L°°(0,T; V). 
Then, we obtain the analogue for p , which implies f 1 1.11) and (II. 3p 2 • ■ 

3 Boundedness, uniqueness, and stability 

In this section, we derive results concerning boundedness, uniqueness and stability of 
the solutions to system fll.ip - fll.4p . With respect to boundedness, we have the following 
result, which resembles Theorem 2.3 in [5]. 

Theorem 3.1 Suppose that (A1)-(A4) are fulfilled, and suppose that the following 
conditions are satisfied: 

(A5) p G £°°(fi), inf p (x) > 0, sup p (x) < 1. 

xen x€n 

(A6) mgL°°(S). 

Then any solution (p, p) of f ll.ip -f lT~4")) fulfilling fl2.3p - fl2.6p also satisfies 

H < /A P > P* 1 and p < p* a.e. in Q (3.1) 

for some constants p* > and p* ,p* G (0, 1) that depend on the structure of the system 
and T , on the initial data, and on an upper bound for the L°° norm of u , only. 

Proof. Let us just show the boundedness of p and the first estimate (13. lft : the results for 
p then follow in exactly the same manner as in the proof of Theorem 2.3 in [5]. Also the 
result for p follows - up to some changes that are necessary due to the different boundary 
condition for p - by the same chain of arguments as in the proof of Theorem 2.3 in [5]; but 
since this proof does not seem to be standard, we provide it for the reader's convenience. 
So let (p,p) be any solution to the system (fT!ft - (fL"4"]L (J22])-(|OI]). We set 

$0 := max{l, ||^ ||L°°(n) , IM|l°°(e)} > 

choose any k G 1R such that k > $0; and introduce the auxiliary function \k £ L°°(Q) 
by putting, for almost every (x, t) G Q , 

Xk(x, t) = 1 if p(x,t) > k, and Xk(x,t) = otherwise. 
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Then, we test (II. ip by (p — k) + . We obtain, for any t G [0, T] , 



c 



+ |2 



|V(/i- k) + \ 2 dxds 



JQ 



+ a (p — u) (p — k) + da ds 



o Jr 

p t — k) + \ 2 dxds — / / p t p (p — k) + dx ds 
o ./n Jo in 

= —k / p t (p — k) + dx ds . 
Jo Jq 

Now observe that a and p are nonnegative and that, by definition of k, 

a (p — u) (/i — k) + = a (\(p — k) + \ 2 + (k — u) (p — k) + ) > a.e. in Q. 
Hence, using Holder's inequality, we obtain from the above equality the estimate 



dl(M(*)-*0 + l& 



|V(/i- k) + \ 2 dxds 



o Jn 



< k / ||Xfc(s)|| L 7/2 (n) ||pt(s)|| L i4/3 (n) ||(//- /c) + (s)|| L2(f1) (is 

Jo 

whence, using the Gronwall-Bellman lemma as in [21 Lemma A. 4, p. 156], 



£|IGu-^) + llco ( [o,T];H) + / / |V(p-A;) + | 2 dx^ V 

Jo Jn 

fc f T 

" ~/| II Pt II L 7 /3(o,T;L 14 /3(C)) IIXfc||L 7 /4(0,T;L 7 /2(n)) • 



(3.2) 



Next, we apply the continuity of the embedding V C L 6 (f2) and the interpolation in- 
equality ( II. 6p with p — 2, q = 6, r = 14/3, and 9 = 1/7. It follows that 



rT 1/3 2 

|Pi||L 7 /3(0,T;L 14 /3(Q)) < (/ || Pt if) || i2(fi) || p t (t) || L 6 (f2) ( 



< II ll 1/7 

— IIP*IIl°°(0,T;#) 



\ 3/7 

IM*)II dt <C\\p t 



1 6/7 



Il 2 (o,t ; v) — ' 



where -Do is a positive constant depending only on the data of the problem. Moreover, 
we have 



|Xfc||L 7 / 4 (0,T;L 7 / 2 (n)) 



o x Jn 



T 
1/2 



\ X k{x,t)\ 7/2 dx) dt 



1/2 



4/7 



I II 8 / 7 

l^ fc llL 2 (0,T;L 4 (Q)) 
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Hence, we can infer from ( 13 .2p that for every k > $o it holds the inequality 

110- k) + \\ < kD x ||Xfc|lK 7 (0 ,T;L4(^)) » (3-3) 
where Z?i = Do/ Min{e, 1} , and where the norm ||| ■ ||| is defined by 

lull 2 := max \\v(t)\\ 2 H + [ [\Vv\ 2 dxdt Wv G C°([0, Tl; if) n L 2 (0, T; V) . 
tepVTI io in 

Moreover, owing to the continuity of the embedding V C L 4 (f2) , there is some constant 
D 2 > , which only depends on f2 and on T , such that 

\Mmo,T;LHn)) <D 2 \lv\\ *iveC\[^T]-H)r\L 2 ^T-V). (3.4) 

At this point, we select a strictly increasing sequence {kj} depending on a real parameter 
m > 1 as follows: 

kj := Af (2 - 2~- J ) for j = 0, 1, . . . , with M := m % . (3.5) 

Note that ko = M > $ and linij_ ) . 00 fcj = 2M . Then, owing to (13. 3p and (13.41) . it is not 
difficult to check that 



(k j+ i - kj) ||Xfc J+ ilU 2 (0,T;L4(n)) < || CM — ^') + IU 2 (0,T;L4(n)) 

1 1 8/7 
k i\\L 2 (Q,T;L A {Q)) 



< D 2 - kj) + \\ < kj D 1 D 2 \\x kj Wl-L T :LHa)Y ( 3 - 6 ) 



Therefore, if we set 

S j := \\Xkj ||L2(0,T;Z,4(n)) for j = 0, 1, . . . , 

then we have 

S j+1 < kj D 1 D 2 sf < AD.D^sf for j = 0, 1, . . . . 

Using [12, Lemma 5.6, p. 95], we can conclude that Sj — > as j — > oo, provided that 

S = \\Xk \\mo,T;LHn)) < (4D lJ D 2 )- 7 2 - 49 . (3.7) 
Now recall that Xk = Xm and, owing to (]3.5j) . M > $ and m = M/$ . Also, 

Ym = 1 < if a > M, and Ym = otherwise. 

M - $ 

Therefore, using (I3.3P and ( 13. 4 p with = ko = M , we find that 

*o< ^^^IIO-^IUw;^)) < M^I^- $ o) + lll 
. DiD 2 ,, ||8 /7 ^ DiD 2 i.s is 

— m _ I IIA<JJt)ll -L J (0,T;L 4 (f2)) — m _ ^ I 

We are now in a position to choose m := 1 + D 1 D 2 |fi| 2/7 T 4/7 (4 £>i D 2 ) 7 2 49 . Then, m > 1 
and (13. 7p is satisfied. Consequently, 

||X2Af||L a (o,T;£*(n)) = lim = 0, 
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due to Beppo Levi's Monotone Convergence Theorem. This implies that p < 2 M a.e. in 
Q , and the boundedness of p is proved. B 



Now that the boundedness condition (13.1 p is shown, we can prove the following uniqueness 
and stability result, which corresponds to Theorem 2.2 in [5]. 

Theorem 3.2 

(i) Suppose that (Al)-(A6) are fulfilled. Then the system fll.ip - fll.4p has a unique solution 
(p, p) satisfying fl2.3l) - fl2.6p . 

(ii) Suppose that (Al) ; (A2) ; (A4) and (A5) are fulfilled and that the functions U\, u 2 
satisfy the conditions (A3) and (A6). Moreover, let (pj,pj) be the solutions to fll.ll) - fll.4p 
corresponding to Ui , i = 1, 2 , and u := U\ — u 2 , p := P\ — P2 and p := p\ — p 2 . Then 
we have, for every t G [0, T] , 



max 

0<s<t 



(\\Ks)\\h+\\p(s)\\V) + / / (\Ms)\\ 2 v + \\ P t(s)\\ 2 H + \\p(s)\\ 2 w ) ds 



Jn 



< K{ I \\u(s)\\l 2(r) ds, (3. 



n 



with a constant K\ > that only depends on the data of the system. 

Proof. Obviously, the assertion (i) follows directly from (ii). So we only need to 
show (ii). To this end, observe that by Theorem 3.1 there are constants M > and 
< r* < r* < 1 such that < p,i < M and r* < pi < r* a.e. in Q , for % — 1, 2 . 
Moreover, the function r 4 r - /'(?")> r * — r — r * 1 nas a Lipschitz constant L > 0. 
Next, we observe that the pair (p, p) is a solution to the system 

(e + 2p 1 )p t + 1pp2,t + (J>Pi,t + H2Pt ~ Ap = a.e. in Q, (3.9) 
Spt - Ap = p - U\pi) - f{P2)) a.e. in Q, (3.10) 

— — = , — — = a(u — p) a. e. on E, (3.11) 
on on 

p(x,0) = p(x,0) = 0, for a.e. xGll. (3.12) 

Now observe that 2p x pp t = (pip 2 ) 4 — p 2 P\,t- Hence, if we test ( 13. 9 p by p then we 
obtain, using Young's inequality, that for every t G [0, T] it holds 

[ (| + pi(t)) p 2 {t)dx + / / |Vp| 2 dxds + / (a\p\ 2 dads 

<cf j\u\ 2 dads+ [ j^-\p\ 2 dads 
Jo Jr Jo Jr 2 



Jn 



|p| (2 |p| |p 2)t | + \p 2 \\pt\)dxds . (3.13) 



We have, owing to the continuity of the embedding if 1 (f2) C L 4 (f2) and to Young's 
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inequality, 

2 \fi\ \p\ \p, 2 ,t\ dxds <C \\fi2,t( 8 )\\H IK s )IU 4 (fi)l|p( s )IU 4 (n) ds 




o Jn 



< 7 / \\p(s)\\ 2 v ds + - f \\p 2 , t (s)\\ 2 H \\p(s)\\ 2 v ds, (3.14) 
Jo 7 Jo 

where, owing to ( 12. 4p . the mapping s h- > \\p>2,t(s)\\% belongs to L 1 (0,T). Moreover, we 
also have fi 2 £ L°°(Q) , and thus 

\p\ \p 2 \ \pt\ dxds <C \\pt(s)\\H\\ld(s)\\ H ds 
'oJn Jo 

< 7 [ t \\pt(s)\\ 2 H ds + - fwn^Wlds. (3.15) 
Jo 7 Jo 

Next, we add p on both sides of Eq. (I3.10P and test the resulting equation by p t . Invoking 
Young's inequality, it is easily seen that, for every t e [0,T] , 



t 

2 J„ , II „U\\\2 



$ \\pt(s)\\ 2 H ds + \\p(t)\\ 2 v 
Jo 

< 7 fht{s)f H ds + - j\\\p{s)\\ 2 H + L 2 ||p( S )||i)^. (3.16) 
Jo 7 Jo 

Now we can combine (I3.13I) - (I3.16I) . Choosing 7 > sufficiently small, and applying 
Gronwall's lemma, we see that (13. 8p is satisfied. B 



The stability estimate ( 13. 8 p can be improved if further regularity is assumed for / . The 
following result is a counterpart of Lemma 3.1 in |7]. We remark at this place that ( 12. 3p 
implies, in particular, that p is weakly continuous as a mapping from [0,T] into W , 
which justifies the formulation of the estimate (3.17) below. 

Theorem 3.3 Suppose that the assumptions of Theorem 3.2, (ii) are satisfied, and 
assume that 

(A7) /GC 3 (0,1). 

Then we have, for every t £ [0,T] , 

(IK* +IM*)|| 2 y + Wp(s)Ww) 



max 

0<s<t 



+ / (\\Pt(s)\\ 2 H + \\pt(s)\\ 2 w ) ds 
Jo 

< K*2 {lK0)||i 2(r) + / (\Hs)\\ 2 L2(r) + \\u t (s)\\ 2 L2{r ^ ds\ (3.17) 



with a constant K\ > that only depends on the data of the system. 

Remark 3.4 We note that ||w(0)||| 2 (r) ^ max{l,t} where I(t) denotes the last 
integral of 03.17p . It follows that ||^(0)|| 2 2 ( r ) can be dropped if one pretends (13.171) just 
for t = T. 
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Proof of Theorem 3.3. We closely follow the lines of the proof of Lemma 3.1 in [JJ. Since 
the proof given there carries over to our situation with minor changes, we can afford to be 
brief. First, observe that by Theorem 3.1 there are constants M > and < r* < r* < 1 
such that < /ij < M and r* < pi < r* a. e. in Q, for i = 1,2. Next, we recall that 
the pair (p, /i) is a solution to the system fl3.9p - fl3.12p . We test Eq. (13. 9p formally by ji t . 
It then follows, with the use of Young's inequality, that 



f \\fi t (s)\\ 2 H ds + l\\Vfi(t)\\ 2 H + ^ j^a\p(t)\ 2 da 



< / a u pt da ds 
Jo Jr 

+ / / (2\p\ \ii 2 ,t\ + M |pi,t| + lA*a| \fh\) \nt\dxds. (3.18) 
Jo Jn 

Now, by virtue of integration by parts with respect to t, and invoking (13.81) . Young's 
inequality and the trace theorem, 

t p 

a u fi t da ds 

o Jr 

< 7 / a\fi{t)\ 2 da + — f |^(t)| 2 da + — [ f \u t \ \fi\ dads 
Jr 1 Jr 7 JoJr 

<C 1 \\fi{t)\\ 2 v + - f \u(0)\ 2 da + - f f(\u\ 2 + \u t \ 2 )dads. (3.19) 
1 Jr 1 Jo Jr 

Employing almost exactly the same arguments as in the proof of Lemma 3.1 in [7] (the 
minor necessary changes are left as an easy exercise to the reader), and taking advantage 
of (13. 8p . we conclude the estimate (where 7 > is arbitrary) 

t f 

(2|p| \fi 2 ,t\ + M \pi,t\ + I ^2 1 \Pt\) \fk\ dxds 




'0 JO. 

< 37 f \\pt(s)\\ 2 H ds + - flMsWvW^fyds 
Jo 7 Jo 

+ - ril^,t(a)lllrl|Ap(s)||lr^ + - / / \u\ 2 da ds . (3.20) 
1 Jo 1 Jo Jr 

Next, we test f)3.10p formally by —Ap t . By the same token as in the proof of Lemma 3.1 
in [7J, we deduce for arbitrary 7 > the estimate 



6 J q \\Vp t (s)\\ 2 H ds + i ||M*)lli < 1 J q W^Wnds 

+ - f (1 + ||p 2 , t (s)|| 2 y ) \\Ap(s)\\ 2 H ds + C [ [ \u\ 2 dads. (3.21) 
7 Jo Jo Jr 

Now observe that, owing to Theorem 2.1, the mappings s ||pi,t(s)|| 2 / , i — 1,2, and 
s ^ llA i 2,t('5)|| 2 ^ all belong to L 1 (0,T). Hence, combining the estimates fl3.18p - fl3.2H) . 
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adjusting 7 > sufficiently small, and invoking Gronwall's lemma, we can conclude that 
for every t G [0,T] it holds 



f (l|Vpt( S )||| + \\m{s)f H )ds + max (\\n{s)\\l + ||p( a )||^) 

JO 0<s<t 

< C |||u(0)||i 2(r) + J J (\u\ 2 + \u t \ 2 ) dads^ . (3.22) 

Next, we formally differentiate (13.101) with respect to t, and obtain 

5p tt - Ap t = nt- /" (pi) Pt - (f"( P i) - f" (pa)) P2,t , (3.23) 



with zero initial and Neumann boundary conditions for p t . Hence, testing (I3.23P by p t 
invoking Young's inequality, and recalling (13.81) and (13.221) . we find that 



5 \\Pt(t)\\ 2 H + JjV Pt (s)\\ 2 H ds 

< C \\\u(0)\\ 2 L , {r) + I [ (\u\ 2 + \u t \ 2 ) dads 



Jr 



+ ff \p2,t\\J"{pi)-f"{p2)\\pt\dxds. (3.24) 
Jo Jn 



Moreover, using Holder's and Young's inequalities, (A7) and (13. 8p . we see that 



/ 

\p2,t\ \f"(pi) ~ f"(p2)\ \pt\ dxds 

'0 Jn 




< C / \\p2,t{s)\\ L A m \\p(s)\\ L i {n )\\pt(s)\\ H ds 

Jo 

< c(J*\\p t (s)\\ 2 H ds + m^J\p(s)\\ 2 v j\\p 2tt (s)\\ 2 v ds) 

< C [ [ {\u\ 2 + \u t \ 2 )dads. (3.25) 

Jo Jr 



Finally, we test (I3.23P by —Ap t . Using Young's inequality and (I3.22p . we find that 
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{ \\V Pt (t)\\ 2 H + J*\\A Pt (s)\\ 2 H ds 

< 7 f \\Ap t (s)\\ 2 H ds + -{|K0)||£ 2(r) + ff {\u\ 2 + h| 2 ) dad 

Jo 7 I Jo Jt 

+ / / \p2,\\f"(pl)-f"(p2)\\^Pt\dxds 

Jo Jn 

< 2 1 j\\Ap t {s)f H ds + ^{||n(0)||i 2(r) + ^(|n| 2 + \u t \ 2 ) dads} 

C f l 
+ — max ||p(s)||^ / \\p 2 ,t{s)\\yds 
7 °<*<* Jo 

< 2 7 f \\Ap t (s)\\ 2 H ds 

Jo 

+ ^ {lk(0)||| 2( r) + ^ (kl 2 + N 2 ) dtrdsj . (3.26) 

Choosing 7 > appropriately small, we can infer that the estimate (I3.17P is in fact true. 
This concludes the proof. a 

4 An optimal boundary control problem 

In this section, we consider the following optimal boundary control problem: 
(CP) Minimize the cost functional 



J(u,p,n) = ^J \p(x,T) - p T (x)\ 2 dx + y J J 



T 

2 



\u(x, t)\ da dt 



+ y / [ Hx,t)-p T (x,t)\ 2 dxdt (4.1) 



Jn 



subject to the state system fll.ip - fll.4l) and to the control constraints 

u G W a d := {v G ^(O.TjL^r)) nL°°(E) : U x < v < U 2 a. e. on E, 

||^IU 2 (0,T;L2(r)) < R] • (4.2) 

In this connection, we require that the hypotheses (Al)-(A7) be satisfied. In addition, 
we postulate: 

(A8) R > 0, ^ > 0, % = 1,2, p T G L 2 (Q), p T G L 2 {Q), U U U 2 G 
L°°(E) , and there are constants < < u* < +00 such that 

< U\ < U 2 < u* a. e. on E . (4.3) 
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In what follows, we denote 

X := H\0, T; L 2 (T)) n L°°(£) , \\u\\ x := |Mk(o,T;£ 2 (r)) + IMU«(E) , 

where || • ||_f/i(o,T ; L 2 (r)) denotes the standard norm in if 1 (0, T; L 2 (r)) . Obviously, W ac j is 
a nonempty, bounded, closed and convex subset of X , and U a d is contained in the open 
set IA C X given by 

f 1 3 

U:=<vEX: - < ess inf v , ess sup W < - w* , ||ut||L 2 (o,T;L 2 (r)) < -R+l 

By Theorem 3.3, the control-to-state mapping u S(u) := (p, p) is Lipschitz continuous 
as a mapping from the set U G X into the space 

(^(0, T; W) n ^([0, T]; V)) x (/^(O, T; if) H C°([0, T]; V)) . 

We may without loss of generality assume (by possibly taking a larger ) that (13. 17)) is 
valid on the whole set U with the same constant K\ > . It also follows from Theorem 3.1 
that there exist constants p* > and < r* < r* < 1 such that for every u G U it holds 

< p < p* and < r* < p < r* < 1 a. e. in Q, (4.4) 

where (p, p) = S(u) . Moreover, a closer inspection of the proof of Theorem 2.1 reveals 
that there is a constant > such that we have, for any wGW, 

\\p\\w 1 >™(o,T;H)nH 1 (o,T;V)nL°°(p,T-,w) 

+ H^II J H' 1 (o,T;H)nc ([o,r];y)nL 2 (o,T;H 3 / 2 (Q))nL° o (Q) < • (4-5) 

Remark 4.1 Thanks to (JOJ) and to / G C 3 (0,1), it holds /'(p) G Also, 
by the embedding V C L 6 (fi), we have p G C°([0, T]; L 6 (fi)) . Notice also that (IZ3JI 
implies, in particular, that p is continuous from [0,T] to H S (Q) for all s < 2; thus, 
since H S (Q) C (7(1]) for s > 3/2, we also have p G C(Q). Therefore, possibly choosing 
a larger constant K£ , we may without loss of generality assume that 

IIpIIc(Q) + IM|c°([0,T];L6(n)) + ||pt|U 2 (0,T;Le(n)) < K3 VllGW. (4.6) 

Remark 4.2 The mathematical literature on control problems for phase field systems 
is scarce and usually restricted to the so-called Caginalp model of phase transitions (see, 
e.g., [UJ, [9], [10], [18], and the references given there). More general, thermodynamically 
consistent phase field models were the subject of p2j. In [7], the present authors investi- 
gated a control problem for the system (ll.ip - (]1.4j) with distributed controls. Since many 
of the arguments employed in [7J carry over to the boundary control considered here, we 
can afford to be sketchy in some of the proofs in the following exposition. 
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4.1 Existence 

We begin our discussion of the control problem (CP) with the following existence result: 

Theorem 4.3 Suppose that the conditions (Al)-(A8) are satisfied. Then the optimal 
control problem (CP) has a solution u £ 

Proof. Let {u n } C W a( j be a minimizing sequence for (CP), and let {(p n , p n )} be the 
sequence of the associated solutions to fll.ip - fll.4l) . We then can infer from (14.51) the 
existence of a triple (u,p,p,) such that, for a suitable subsequence again indexed by n, 
we have 

u n — y u weakly star in H\0, T; L 2 (T)) n L°°(X)), 

p„ ^p weakly star in W 1,oo (0, T; H) D H l {Q, T; V) D L°°(0, T; W), 

/i n ^p weakly star in H 1 (0, T; if) fl L°°([0, T];V) C\ L 2 (0, T; H 3 / 2 (£l)). 

Clearly, we have that u £ W a d . Moreover, by virtue of the Aubin-Lions lemma (cf. [T4"l 
Thm. 5.1, p. 58]) and similar compactness results (cf. [XTf Sect. 8, Cor. 4]), we also have 
the strong convergences 

p n -»■ p strongly in C°([0, T]; if s (fi)) for all s < 2, 
strongly in C°([0, T); H) D L 2 (0, T; V). 

From this we infer, possibly selecting another subsequence again indexed by n, that 
Pn P pointwise a. e. (actually, uniformly) in Q . In particular, r* < p < r* a. e. in 
Q and, since / £ C 2 (0, 1), also f'(p n ) f'(p) strongly in L 2 (Q). Now notice that the 
above convergences imply, in particular, that 

p n ^p strongly in C°([0,T];L 6 (fi)), 
d tPn -> d t p weakly in L 2 (0, T; L 4 (fi)), 
yU„ — > p strongly in L 2 (0, T; L 4 (fi)), 
dtp n — > d t p weakly in L 2 (Q). 

From this, it is easily verified that 

Pn d t Pn -> P d t p weakly in L 1 (0, T; if) , 
pndtPn^pdtp weakly in L 2 (0, T; L 3/2 (Q)). 

In summary, if we pass to the limit as n — > oo in the state equations fll.ip - fll.4l) written for 
the triple (u n , p n , p n ) , we find that (p, p) = S(u) , that is, the triple (u, p, p) is admissible 
for the control problem (CP). From the weak sequential lower semicontinuity of the cost 
functional J it finally follows that u, together with (p,p) = S(u) , is a solution to (CP). 
This concludes the proof. a 

Remark 4.4 It can be shown that this existence result holds for much more general cost 
functionals. All we need is that J enjoy appropriate weak sequential lower semicontinuity 
properties that match the above weak convergences. 
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Remark 4.5 Since the state component p is continuous on Q, the existence result re- 
mains valid if suitable pointwise state constraints for p are added (provided the admissible 
set is not empty). 

4.2 Necessary optimality conditions 

In this section, we derive the first-order necessary conditions of optimality for problem 
(CP). To this end, we first show that the control-to-state operator S : u y-¥ (p, p) is 
Frechet differentiable as a mapping from U C X into the Banach space (y, \\ ■ \\y) , 
where 

y := (H 1 (0,T;H)nC ([0,T};V)f]L 2 (0 J T;W)) 
x (C°([0,T];H)nL 2 (0,T;V)). 

4.2.1 The linearized system 

Let u G U and h G X be given and (p, p) = S(u) . As a preparatory step, we consider 



the following system, which is obtained by linearizing the system fll.ip -f TT~4"|) at (p, p) : 

(e + 2p)r] t -Ar] + 2p t ^ + p^ + p t 7] = a.e. in Q, (4.7) 

5£ t -A£ = -f(p)£ + 77 a.e. in Q, (4.8) 
<9£ drj 

— — = , — — = a(h — n) a. e. on E, (4.9) 
on on 

£(x, 0) = rj(x t 0) = for a.e. i£fl. (4.10) 

We expect for the Frechet derivative DS(u) at u (if it exists) that (£, r/) = DS(u)h, 
provided that (I4.7p - fl4.10p admits a unique solution (£,?/). In view of (12.31) . (12.41) . and 
(13.11) . we can guess the regularity of £ and rj : 

£ G ^(0, T; if) n C°([0, T]; V) D L 2 (0, T; W) H f °°(Q), (4.11) 
G ^(O.T; If) nC°([0,T];F) nf 2 (0,T; Jf 3/2 (fi)). (4.12) 



Notice that also in this case we cannot expect that p(t) G H 2 (Q) a.e. in (0, T) due to the 
low space regularity of h , and we could repeat Remark 2.2 here. Nevertheless, if (14. lip and 
(I4.12p hold, then the collection of source terms in (14.71) . i. e., the part — 2 p t £ — p £ t — p t rj , 
belongs to L 2 (Q), whereas the regularity (I4.12p for p allows us to conclude from (14. 8p 
that also £ G C(Q) (by applying maximal parabolic regularity theory, see, e.g., [H Thm. 
6.8] or [El Lemma 7.12]). 

In fact, £ is even more regular: indeed, we may differentiate (14. 8p with respect to t to 
find that 

6{ u - A& = -f'"(p) p t { - f (p) t t + rk , (4.13) 

with zero initial and Neumann boundary conditions for £ t . Since the right-hand side 
of (I4.13P belongs to L 2 (Q), we may test by any of the functions £ t , £ tt , and — A£ t , to 
obtain that even 



£ G Jf 2 (0,T; i/)nC 1 ([0,T];y)nff 1 (0 1 T;^) 



(4.14) 
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Notice, however, that this fact has no bearing on the regularity of rj , since the coefficient 
fit in (14 .7p only belongs to L 2 (Q) . 

The following well-posedness result resembles Proposition 3.2 in [7]. 

Proposition 4.6 Suppose that (A1)-(A8) are fulfilled. Then the system fj4T7|) - fj4TT0|) 
has a unique solution (£,??) satisfying (14.121) . (14. 14j) . and 

H^I|H 2 (0,T; J H')nC 1 ([0,T];y)n J H" 1 (0,T;Ty) + H ? ?l| J H' 1 (0,T; J H')nC ([0,T];y)nL 2 (0,T; J H' 3 / 2 (n)) 

< Kl\\h\\ W]LHr)) , (4.15) 

with a constant K\ > that is independent of the choice of u £ W and h £ X . 

Remark 4.7 It follows from Proposition 4.6, in particular, that the linear mapping 
h (£,rj) is continuous from X into y. 

Proof. We follow the lines of the proof of our previous existence results and proceed in a 
series of steps. 

Step 1: Approximation. As in the proof of Theorem 2.1, we use an approximation 
technique based on a delay in the right-hand side of (]4.8j) . To this end, for r > we 
resume the definition of the translation operator T T : L 1 (0,T;if) — > L 1 (0,T;if) by 
putting, for every v £ L 1 (0,T; H) and almost every t £ (0, T) , 

(T T v)(t) = v{t -r) if t > r, and (7»(t) = if t < r. (4.16) 

Notice that, for any v £ L 2 (Q) and any r > 0, we obviously have \\%v \\l 2 (Q) < IM|l 2 (Q) ■ 
Then, for any fixed r > , we look for functions (£ r , r] T ) , which satisfy (14.111) and (14. 12j) 
and the system: 

{e + 2p)r]l -Ar] T + 2fi t C + + PtV T = a.e. in Q, (4.17) 
5Q - AC + f'(p)C = T rV T a.e. in Q, (4.18) 

^1 = 0, ^L = a (h-ri T ) a.e. on S, (4.19) 
an on 

C(x, 0) = r] T (x, 0) = for a. e. x £ Q. (4.20) 

Precisely, we choose for r > the discrete values r = T/N , where £ IN is arbitrary, 
and put t n = nr , < n < N , and J n = (0,t n ) . For 1 < n < N , we solve the problem 

(e + 2p)r] n j - Ar] n + 2fi t £ n + fi£n,t + PtVn = a.e. in fix 7 n , (4.21) 

-Tp 1 = a(h — rj n ) a. e. on T x J n , rj n (x, 0) = for a. e. x £ f2, (4.22) 

5 Cn,t - A£ n + f"(p) C = T T rin a. e. in Q x / n , (4.23) 

^P = a.e. on Tx 4, ^ n (x,0) = for a.e. x £ fi, (4.24) 
an 

where the variables r\ n and ^ n , defined on , have obvious meaning. Here, T T acts on 
functions that are not defined on the entire interval (0, T) ; however, for n > 1 it is still 
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defined by (I4.16P , while for n = 1 we simply put T T ri n = . Notice that whenever the 
pairs (Ck,Vk) with 

£ fe G H\l k ; H) D C°(/ fc ; V) PI L 2 (I k ; W) PI C(lT^4), (4.25) 

% e ^(4; if) n C°(4; VO n L 2 (4; # 3/2 (^)), (4.26) 

have been constructed for 1 < k < n < N , then we look for the pair (^+1,77^+1) that 
coincides with {i n ,rj n ) in J n , and note that the linear parabolic problem (I4.23p . (I4.24p 
has a unique solution f n+1 on Q x I n+1 that satisfies (I4.25P for k = n + 1 . Inserting f n+1 
in (14.211) (where n is replaced by n + 1 ) , we then find that the linear parabolic problem 
f)4.2ip . (I4.22p admits a unique solution rj n+ i that fulfills (14.261) for k — n + 1 . Hence, we 
conclude that (f r ,*7 T ) = (6v,w) satisfies fl4TT7]) - flCT]) . and fICTD . 

Step 2: A priori estimates. We now prove a series of a priori estimates for the functions 
(£ T ) r / r ) • I n the following, we denote by Ci (i £ IN) some generic positive constants, which 
may depend on e, 5, p*, p*, p*, T, K*, K%, , but not on r (i.e., not on N). For the 
sake of simplicity, we omit the superscript r and simply write (£,77) . 

First a priori estimate. Observe that 2pr]r]t = {prj 2 ) t — p t if. Hence, testing (I4.17P by 
t], and invoking (14.191) and Young's inequality, we have, for < t < T , 

J (| + P(t)) V(t? dx + J \\Vr](s)\\ 2 H ds + J^J arfdads 

< I I ' -rf dads + C x [ f h 2 da ds 
Jo Jr 2 Jo Jv 

+ 2 \P*\ If I \v\ dxds + / / \p\ \£ t \ \rj\ dxds. (4.27) 
Jo Jn Jo Jn 



For any 7 > 0, we have, by Young's inequality and (14. 4p . that 

|/2| |f t | \r)\dxds < \\fl\\ L °°(Q) / \\v(s)\\HUt(s)\\ H ds 



Jn 



\Us)\\ 2 H ds + °* / \\v(s)\\ 2 H ds. (4.28) 
7 Jo 



Moreover, 



\p t \ |f I \v\dxds < / \\pt(s)\\ H U(s)\\ L 4 in) \\rj(s)\\ L 4 {n) ds 



Jn 



2 j„ , ^3 



< 7 / hOOII^ + ^ / ||/2 t ( S )||^ \\m\\ 2 vds. (4.29) 
Jo 7 Jo 

Notice that, by virtue of (14. 5p . the mapping s 1— >■ ||/it(s)||# is bounded by a function in 
^(0,T). 

Next, we add f on both sides of Eq. (14.181) and test the resulting equation by f t . On 
using Young's inequality again, we obtain: 

lJ*\Ms)\\ 2 H ds + ~ (||^)||| + ||V£(t)||!r) 

< c 4 f r 11^)11^^+ r 11^)11^^). (4.30) 
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Combining the inequalities (I4.27p -f l4.30p . and choosing 7 > sufficiently small, we con- 
clude from Gronwall's lemma that 







{\Ut)\\H + Ht)\\v)dt + max (||£(i)f v + \\ v (t)\\%) 

<C 5 [ [ \h\ 2 dadt. (4.31] 
Jo Jv 



Thanks to (14.191) , we may also infer (possibly by choosing a larger C5 ) that 



U(t)\\w < C 5 {\\Aat)\\ 2 H + ^J r \h\ 2 dadt) for all te [0,T]. (4.32) 
Second a priori estimate. We test (14.171) by r] t and apply Young's inequality in order 



to obtain 



ej \\Vt(s)\\ 2 H ds + ^\\Vv(t)\\ 2 H + J\ \ri(t)\ 2 da < J Jahrftdads 

(2 \p t \ |f | + |/2| |6l + \pt\ \v\ ) \Vt\ dx ds. (4.33) 



'0 Jn 

By (14. 4p . we can infer from Young's inequality that 



f [ \fi\\tt\\Vt\dxds < 7 fht{s)\\ 2 H ds + ^ f \\Us)\\ 2 H ds. (4.34) 
Jo Jn Jo 7 Jo 



Moreover, by virtue of Holder's and Young's inequalities, 



\pt\ \v\ \Vt\ dx ds 



Jn 

2 > , °7 



' C 7 " 



<1 I \\Vt(s)\\ H ds + — \\pt(s)\\ LHn) \\v(s)\\ LHQ) ds 
Jo 7 Jo 



t n rt 



<7 / \\Vt(s)\\ 2 H ds + — / \\pt(s)\\vMs)\\ 2 v ds. (4-35) 
Jo 7 Jo 

Observe that by (14. 5 p the mapping s ||pt(s)||y is bounded by a function in L 1 (0,T) . 
Also, we have, owing to the continuity of the embedding W C L°°(Q) and (14.321) . 



/ / 2 \n t \ |£| \rj t \ dxds 
Jo Jn 



< 7 I \\Vt(s)\\ 2 H ds + ^£ f UMf H U{s)\\l^n)ds 
Jo 7 Jo 



< 1 fht{s)\\ 2 H ds + ^^^\h\ 2 dadt 

+ j)\Ms)\\ 2 H \ms)\\ 2 H ds), (4.36) 
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where, owing to (14. 5p . the mapping s H- ||/it(s)||# is bounded by a function in L 1 (0,T) . 
Finally, we employ integration by parts, Young's inequality, (I4.3ip . and the trace theorem 
to obtain 



ahrjtda ds 



o Jr 



< / a | h(t) | \r](t) | da + / / a \h t \ \r)\ dads 
Jr Jo Jr 

< J^\v(t)\ 2 da + Cn\\h\\ 2 H1 ^ T . L2{r)) . 
Next, we formally test (I4.18P by — A^ t to obtain, for every t G [0, T] , 

s[ t \\V^)\\ 2 Hds + h\Am\ 2 H= f! (-(TrV)-f"(p)0 Htdxds. 

Jo 1 Jo Jn 

Now, by virtue of (I4.3ip and invoking Young's inequality, we have 

n(T T r]) A£ t dx ds 

< [\(%v)(t)\ \A£(t)\dx+ I I \d t {T T r))\ \A^\dxds 
Jn Jo Jn 

< l\\m)\\H + C 12 £ £\h\ 2 da dt 

+ 7j\\Vt(s)\\%ds+^j]\Am\\ 2 H ds. 



Moreover, it turns out that 



/ / f"(p)ZA£ t dxds 
Jo Jn 



< i\f"(p(t))\m)\mt)\dx 

In 



I I \f"\p)M + f"{p)tt\ \Ai\dxds. 
Jo Jn 

We have, owing to P~4l) and flOB . 

)f"(p(t))\\m\\m)\dx < lm(t)\\ 2 H + c l3 £^\h\ 2 dadt. 



(4.37) 



(4.38) 



(4.39) 



(4.40) 



(4.41) 



Also the second integral on the right-hand side of (14.401) is bounded, since (I4.4p . (14. 5p . 
and (14.311) imply that 



ff \f"\p)Pti + f\p)it\ \Ad\dxds 
Jo Jn 



< C u /(||Pt(«)ll| 4 (n)lie( S )||i 4( n) + liet( S )||lr)^+ /WoOII^ 

Jo Jo 

< Ci5(max ||^)||^||p,( S )ir^ S + £||6( S )||^d S ) +j)ms)\\ 2 H ds 

< Cie I* f \h\ 2 dadt+ []\AC(s)\\ 2 H ds, (4.42) 

Jo Jr Jo 
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thanks to the continuity of the embedding V C L 4 (Q) . Thus, combining the estimates 
(I4.33p — (I4.42p . choosing 7 > sufficiently small, and invoking Gronwall's inequality, we 
can infer that 

[\\\Vt(t)\\ 2 H + \Ut)\\ 2 v)dt + max (||^)l| 2 y + U(t)\\ 2 W ) 

< Cir IWlHi(o,T;£2(r)) ■ (4.43) 

Next, we compare terms in ( 14.1 7\ and, arguing as in the derivation of (I4.33p -( l4\37p . we 
readily find that 

/ \\A v (t)\\ 2 H dt < C 18 \\h\\ 
Jo 

Thus, by owing to elliptic regularity (cf. (I4.19P and Remark 2.2), we conclude that 

/ \\v{t)\\ 2 H 3/2 ia) dt < C 19 ||/i||| 1(0)T . i2(r)) . (4.44) 
Jo 

Finally, we differentiate Eq. (I4.18P with respect to t . We obtain: 

SCu-Mt = d t (r T v)-f m (p)ptC-f"(p)Ct a.e. in Q. (4.45) 

From ( 0]) -( l4~6]) . f )4.43[) and (I4.44p . we can infer that we may test (I4.45P by any of the 
functions £t , — A£ t , and £ tt , in order to find that 

[ T (Uu(t)\\ 2 H + l|A6(*)||^) dt + max \Ut)\\ 2 v < C 20 \\h\\ 2 HH0>T . L2(r)) . (4.46) 

Step 3: Passage to the limit. Let (^ T ,r] T ) denote the solution to the system (14.171) - 
(I4.20p associated with r = T/N , for iV G IN. In Step 2, we have shown that there is 
some C > , which does not depend on r , such that 

ll^ r |lH 2 (o,r;_H")nc* 1 ([o,T] ; y)nH 1 (o,T;iy)nC(Q) 

+ ll 7 7 T |lH 1 (0,T;H)nC ([0,T];1/)nL 2 (o,T;i/3/2( C )) < C || tl \\ H i (o,T;L 2 (r)) • (4.47) 

Hence, there is a subsequence \ such that 

C k ^£ weakly star in H 2 (0, T; H) D W /1,oo (0, T; V) PI i^ 1 (0, T; W), 
7] Tk ^7] weakly star in ^(0, T; H) D L°°(0, T; V) n L 2 (0, T; if 3/2 (fi)). 

(4.48) 

From the trace theorem we can infer that £ satisfies the boundary condition given in 
(14. 9p . while the boundary condition for 77 will be satisfied (either in the variational sense 
or in the sense of the appropriate trace theorem, see Remark 2.2) once we prove that we 
can pass to the limit in the products of ( 14 .7p . as shown below. Moreover, it is easily seen 
that also ( 14. lUp is fulfilled. By compact embedding, we also have, in particular, 

C k -> £ strongly in C(Q), rf k ->> r] strongly in L 2 (Q), (4.49) 
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so that prf t k -> pr] t and p& k both weakly in L 2 (Q) , /"(p)£ Tfc ~> /"(p) £ 

strongly in L 2 (Q) , as well as //t£J fe — > pt^t and pt?7 rfc — > p t 77 , both strongly in L l (Q) . 
Finally, it is easily verified that {T Tk ri Tk } converges strongly in L 2 (Q) to 77 . In conclusion, 
we may pass to the limit as k — >■ 00 in the system (j4.17p - (l4.20p (written for r k ) to find 
that the pair (£,77) is in fact a solution to the linearized system (I4.7p - (l4.10p . 
We now show the uniqueness. If (£1,771), (£2)^2) are two solutions having the above 
properties, then the pair (£, 77) , where £ = £1 — £2 and 77 = 771 — 772 , satisfies (I4.7p - fl4.10p 
with h = 0. We thus may repeat the first a priori estimate in Step 2 to conclude that 
£ = 77 = 0. 

Finally, taking the limit as r \ in (I4.47P and invoking the lower semicontinuity of 
norms, we obtain the inequality f )4.15p . This concludes the proof. 



4.2.2 Frechet differentiability of the control-to-state mapping 

In this section, we prove the following result. 

Proposition 4.8 Suppose that the assumptions (Al)-(A8) are satisfied. Then the 
solution operator S , viewed as a mapping from X to y , is Frechet differentiable on U . 
For any u G U the Frechet derivative DS{u) is for h G X given by DS{u)h = (£,77), 
where (£,77) is the unique solution to the linearized system fl4.7p - fl4.10p . 

Proof. Let u G U be given and (p, p) = S(u) . Since U is an open subset of X , there 
is some A > such that u + h EU whenever h G X satisfies \\h\\ x < A . In the following, 
we consider such perturbations h G X , and we define (p h ,p h ) := S(u + h) and put 

z h := p h - p - 7] h , y h -=p h -p-^ h , (4.50) 

where (£ ft , rf 1 ) denotes the unique solution to the linearized system fl4.7l) - fl4.10l) associated 
with h. Since the linear mapping h 1— > (£ h ,77 /l ) is by Proposition 4.6 continuous from 
X into y, it obviously suffices to show that there is an increasing function g : [0, A] — > 
[0, +00) which satisfies lim,,^ g{r)/r 2 = and 

WSJ ll# 1 (0,T;i?)nC°([0,T];V)nL 2 (0,'7W) \\ Z II C°([0,T];H)nL 2 (0,T;V) 

< 9 (IW|ffi(o,T;i,2(r))) . (4.51) 

Using the state system (ll.2|) - fll.4p and the linearized system (I4.7p - fl4.10p . we easily verify 
that for h G X with \\h\\x < A the pair (y h ,z h ) is a strong solution to the system 

(e + 2p) 4 + p t z h + p y h t + 2p t y h - Az h 

= -2 (p h t - Ji t ) (p h - p) - (pt - p t ) (p h - p) a. e. in Q, (4.52) 
6tf - Ay h + f\p h ) - f{p) - f{p) i h = z\ a. e. in Q, (4.53) 

= 0, = -az h , a.e. on E, (4.54) 

an an 

y h (x,0) = z h (x,0) =0 for a.e. x E tt. (4.55) 
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Notice that 

y h G H 1 ^, T; H) n C°([0, T]; V) D L 2 (0, T; n C(Q), 
^ G ^(0, T; #) n C°([0, T]; V) n L 2 (0, T; ff 3/2 (fi)). 

For the sake of a better readability, in the following estimates we omit the superscript 
h of y h and z h . Also, we denote by Cj (i G IN) certain positive constants that only 
depend on e, 5, p*, p*, p*, T, iT*, ifg, , but not on /i. 

We now add y on both sides of Eq. (14.531) and test the resulting equation by y t . Using 
Young's inequality, we find that for all t G [0, T] it holds 



\z(s)\\ H ds 



^ll^)lll^+^(l|V^)lHr+lb(*)K < \f 
+ Ci f IbOOIIlr^ + C 2 ril(/V)-/ / (p)-r(p)e' l )( S )|^^- (4-56) 



In order to handle the third term on the right-hand side of (I4.56p . we note that the 
stability estimate (I3.17P implies, in particular, that 



\\P ~ p\\l°°(q) < K% IWIh^o.t;!, 2 ^)) > (4-57) 

that is, p h p uniformly on Q as ||/i||,?/i(o,T;L 2 (r)) — ^ 0. Since / G C 3 (0, 1) , we can infer 
from Taylor's theorem and (I4.4p that 

\f{p h )-f{p)-f"{p)e\< max l -^^\p h -p\ 2 +\f"(p)\\y\ onQ. (4.58) 
It then follows from the estimates (13.171) and (I4.56I) - (I4.57I) that 

\[ht{s)\\lds + \\\y{t)\\l < -j\\z{s)\\ 2 H ds + C,j\\y{s)f H ds 

+ C,\\h\\^ T . L , {r)) . (4-59) 

Next, observe that 2pzz t = (p z 2 ) t — p t z 2 . Therefore, testing (I4.52p by z yields for every 
t G [0, T] that 

/ (1 + P^) z2 (t)dx + [ \\Vz(s)\\ 2 H ds + [ [a\z\ 2 dadt 

(p y t + 2 \± t y) z dx ds — 2 / / (pj* — p t ) (p h — p) 2 efcc 
Jn Jo Jn 



./n 



{p\ l - p t ) (fi h - p) z dx ds. (4.60) 



We estimate the terms on the right-hand side of (I4.60p individually. At first, using (I4.4p 
and Young's inequality, we find that 

/ / \m\y t \\z\dxds < 7 / \\y t (s)\\ 2 H ds + ^ / \\z(s)\\ 2 H ds. (4.61) 
Jo Jn Jo 7 Jo 



Colli Gilardi Sprekels 



27 



Moreover, using the continuity of the embedding H l (Q) C L 4 (f2) , as well as Holder's and 
Young's inequalities, we have 



fi t \ \y\ \z\ dxds < 2 / \\pt(s)\\ H \\z(s)\\ L 4(n) \\y(s)\\L*(si)ds 
Jo 

< 7 f \\z{s)\\ 2 v ds + ^ f\\Hs)\\ 2 H \W)f v ds. (4.62) 




o Jn 



7 Jo 

Observe that by (I2.4p the mapping s >->■ ||/it(s)||jj belongs to L 1 (0,T). 

At this point, we can conclude from (13.17P and (I4.57p . invoking Young's inequality, that 

2 \nt — p*t\ \p h — p\ \z\dxds 

'o Jn 

< 2 J - (it)(s)\\ H ||(/-p)(s)|U (n) 

ds + y ||«(s)||^da 

< / \\z(s)\\ 2 H ds + C 8 ||/i||^ ( o,T;L 2 (r)) • (4-63) 
Jo 

Finally, we invoke (I3.17P and Holder's and Young's inequalities, as well as the continuity 
of the embedding if 4 (f2) C L 4 (f2) , to obtain that 

t 

— pt\ \p h — p\ \z\dxds 

'o Jn 

< max \\z(s)\\ H / \\(pt ~ Pt)(s)\\ L4{n) \\(p h - p)(s)\\ LHn) ds 



0<s<t 



n 



C f f 

< 7 max -|k(s)||# + — / ||(p£ - p*)0) ||y ^ / || (// - p)(s) \\ 2 y ds 
u - s -' 7 jo >/o 

- 7 0™ II^Hh + ° W W h WHi{0,T;L*{T)) ■ ( 4 - 64 ) 

Combining the estimates (I4.59I) - (I4.64I) . taking the maximum with respect to t G [0,T] , 
adjusting 7 > appropriately small, and invoking Gronwall's lemma, we arrive at the 
conclusion that (y h , z h ) = (y, z) satisfies the inequality 

\\U \\m{0,T;H)r\C°{[Q,T\;V) ^~ W Z \\C°([0,T];H)nL 2 ([0,T];V) 

< C11 ll^ll^co^CD) ■ (4-65) 
Finally, testing (14.531) by —Ay h , and using (j4.58p . we find that also 

\\V \\l*{0,T;W) - Cl2 IWIffl(0,T;Z,2(r)) • (4.66) 

Therefore, the function g(r) := (Cn + Ci 2 ) r 4 has the requested properties. This con- 
cludes the proof of the assertion. a 
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Corollary 4.9 Let the assumptions (Al)-(A8) be fulfilled, and let u G W a d be an 
optimal control for the problem (CP) with associated state (p,p) = S(u) . Then, for 
every v G U ad , 

fau(v-u)dadt + (p(T) - p T ) £(T) dx + / / fa (p ~ A*r) V dx dt > 0, (4.67) 




o Jr Jn Jo Jn 



where (£, 77) is the unique solution to the linearized system (]4.7p -( l4TT0~]) associated with 
h = v — u . 

Proof. Let v G W a d be arbitrary and h — v — u. Then u + Xh G W a( j for < A < 1 . 
For any such A , we have 

< Jju + Xh, S(u + Xh)) - J(u, S(u)) 
X 

J(u + Xh, S{u + Xh)) - J(u, S(u + Xh)) 



< 



+ 



A 

J(u,S(u + Xh)) - J(u,S(u)) 
X 



It follows immediately from the definition of the cost functional J that the first summand 
on the right-hand side of this inequality converges to f Q f r f3 1 u(v — u) da dt as A \ . 
For the second summand, we obtain from Proposition 4.8 that 

J(u,S(u+ Xh)) - J(u,S(u)) 
hm 

A\0 A 

\p{x, T) — pr(x))£(x, T)dx + / I fa {p - P-t) V dx dt , 



'a Jo Jn 

whence the assertion follows. B 

4.2.3 The optimality system 

Let u G W a d be an optimal control for (CP) with associated state (p, pi) = S(u) . Then, 
for every v G W a d , (I4.67P holds. We now aim to eliminate (£, 77) by introducing the 
adjoint state variables. To this end, we consider the adjoint system: 

- (e + 2p) q t - p t q - Aq = p + fa (p - p T ) a. e. in Q, (4.68) 
dq 

— — = —aq a. e. in S, q(x,T) = for a. e. x G f2, (4.69) 
an 

-Spt -Ap + f"(p) p = pq t -p t q in Q, (4.70) 

^ = on E, 6p(T) = p(T) - p T in Q , (4.71) 
an 

which is a linear backward-in-time parabolic system for the adjoint state variables p and 
?■ 
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It must be expected that the adjoint state variables (p, q) be less regular than the state 
variables (p,p)- Indeed, we only have p(T) G L 2 (Q), and thus (I4.70P and f)4.7ip should 
be interpreted in the ususal weak sense. That is, we look for a vector-valued function 
p G H\0, T; V*) nC°([0, T]]H) flL 2 (0, T; V") that, in addition to the final time condition 
flUEED, satisfies 

(-5p t (t),v) v * tV + f Vp(t)-Vvdx+ [ f"(p(t))p(t)vdx 
Jn Jn 

= I (/*(*) Qt(t) - Pt{t) q{t)) v dx , (4.72) 
Jn 

for every v G V and almost every t G (0,T). Notice that if q G if x (0, T;if) D 
C°([0, T]; V) , then it is easily seen that pq t — p t q G L 3 ^ 2 (Q) , so that the integral on the 
right-hand side of (I4.72p makes sense. On the other hand, if p has the expected regu- 
larity then the solution to (Qgj) . flOT]) should belong to H\Q,T\H) n C°([0,T]; V) n 
L 2 (0,T;tf 2 (fi)). 

The following result is an analogue of Theorem 3.7 in [7]. 



Theorem 4.10 Suppose that u G W a d «s an optimal control for (CP) wit/j associated 
state (p, £t) = iS'(w) . Then the adjoint system fl4.68l) - fl4.71l) has a unique weak solution 
(p,q) withpe H\0,T;V*)nC ({0,T};H)r)L 2 (0,T;V), q G if^O, T; H)f)C ([0, T]; V)f) 
L 2 (0,T; H 2 (Q)) ; moreover, for any v G U & &, we have the inequality 

(3\u{v — u) da dt + / / a g (u — it) da dt > . (4.73) 




Proof. The existence and uniqueness result for the adjoint state variables p and q 
follows using the same line of arguments as in the proof of Proposition 3.6 in [7], with 
only minor and straightforward changes that are due to the different boundary condition 
for q. Now let v G be given. A standard calculation (which can be left as an easy 
exercise to the reader), using the linearized system fl4.7l) - (l4.10l) with h = v — u, repeated 
integration by parts, and the well-known integration by parts formula 

/ {{v t (t),w{t))v* tV + (w t (t),v(t)) v *y)dt= [ (v(T)w(T) - v(0)w(0)) dx 
Jo Jn 

(which holds for all functions v, w G -f^ 1 (0, T; V*) H L 2 (0, T; V) ), yields the identity 
(p(x, T) - pt(x))£(x,T) dx + / f3 2 (p - pr)v dx dt 




n Jo Jn 

T 

' aq(v - u)dadt. (4.74) 
'o Jv 

The variational inequality (14.731) is thus a direct consequence of Corollary 4.9. 
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